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ABSTRACT 

We present IRAM Plateau de Bure Interferometer 3 mm observations of CO 7(1— >0) emis- 
sion in two 24/.tm-selected starburst galaxies in the outskirts (^2-3i? v ; r i a i) of the rich clus- 
ter CI 0024+16 (z — 0.395). The galaxies' inferred far-infrared luminosities place them in 
the luminous infrared galaxy class (LIRGs, Lfir > lO n L ), with star formation rates 
of ^60A/ Q yr _1 . Strong CO 7(1— >0) emission is detected in both galaxies, and we use 
the CO line luminosity to estimate the mass of cold molecular gas, A/(H2). Assuming 
M(R 2 )/L' co = O.8M (Kkm- 1 pc 2 )- 1 , we estimate M(H 2 ) = (5.4-9.1) x 1O 9 M for 
the two galaxies. We estimate the galaxies' dynamical masses from their CO line-widths, 
Mdyn ~ 1-3 x 1O 1O M , implying large cold gas fractions in the galaxies' central regions. At 
their current rates they will complete the assembly of M* ~ 1O 1O A7 and double their stellar 
mass within as little as ^150Myr. If these galaxies are destined to evolve into SOs, then the 
short time-scale for stellar mass assembly implies that their major episode of bulge growth 
occurs while they are still in the cluster outskirts, long before they reach the core regions. 
Subsequent fading of the disc component relative to the stellar bulge after the gas reservoirs 
have been exhausted could complete the transformation of spiral-to-SO. 
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1 INTRODUCTION 

One key requirement for a starburst is the presence of a reservoir of 
dense, cold gas that can be efficiently converted to stars. For galax- 
ies entering rich clusters this is especially important, because they 
are expected to be affected by mechanisms that can remove cold gas 
from the haloes and discs of infalling galaxies (e.g. ram-pressure 
stripping, Gunn & Gott 1972) or prevent further cooling of gas 
within galaxies' dark matter halos (starvation or strangulation, e.g. 
Larson, Tinsley & Caldwell 1980; Bekki, Couch & Shioya 2002). 
This environmental dependence has a profound influence on clus- 
ter galaxies' evolutionary histories, the net effect of which is the 
eventual termination of star formation (there is virtually no resid- 
ual star formation in the cores of local clusters). The observational 
evidence for the gradual truncation of star formation in clusters is 
the conspicuous disappearance of star-forming disc galaxies in the 
cores of rich clusters since z ~ 1. It has been proposed that the 
spiral galaxies in distant clusters must be transforming into passive 
lenticular (SO) galaxies, since the fraction of SOs in rich clusters is 
believed to increase at the same time that the spirals begin to vanish 



(Dressier et al. 1997). However, the detailed nature of the physics 
controlling this evolution is still poorly understood. For example, 
what process is responsible for the transformation of the bulge-to- 
disc ratio of the spiral population? (Kodama & Smail 2001). 

A growing body of observational evidence suggests that dis- 
tant clusters contain significant populations of dust obscured star- 
burst galaxies in their peripheral regions (e.g. Due et al. 2000; Bi- 
viano et al. 2004; Coia et al. 2004; Geach et al. 2006; Marcillac 
et al. 2007; Bai et al. 2007; Dressier et al. 2008; Koyama et al. 
2008; Fadda et al. 2008). Geach et al. (2009) proposed that the 
population of luminous infrared galaxies (LIRGs) residing in the 
'infall' population of z ~ 0.5 clusters could be examples of spi- 
ral galaxies undergoing an episode of bulge growth via circumnu- 
clear starburst, and therefore excellent candidates for the progeni- 
tors of local massive SOs. Geach et al. (2009) estimate that, given 
the cluster LIRGs' current SFRs, 10 10 M Q of stars could be built 
up within a few 100 Myr. Although this simple model suggests that 
these galaxies could evolve onto the bright end of the local cluster 
SO luminosity function by z — 0, we still lack crucial observational 
constraints, such as the likely duration of the starburst, where it oc- 
curs in the cluster, and limits on the final stellar mass of the bulge. 
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In this Letter we refine our model by asking a simple question: 
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Table 1. Details of the two galaxies in our study: MIPS J002621.7+171925.7 and MIPS J002721.0+165947.3 (Geach et al. 2009). We present the results of the 
CO observations: CO(l-O) line widths, luminosities and corresponding H2 gas mass. Mid-infrared observations are from Geach et al. (2006 & 2009), and the 
spectroscopic redshifts are from Czoske et al. (2001) and derived from the [O II] emission. Note the double Gaussian profile of MIPS J00272 1 .0 was a four 
parameter fit, with a single FWHM and amplitude and two velocity offsets. 
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f Offset from [O II] line. MIPS J002721.0 is best fit with a double Gaussian, and we quote the offsets of both peaks. The uncertainty is ~30kms 1 . 
ft Assuming a CO luminosity to total gas mass conversion factor of a = O.SMQfKkm" 1 pc 2 ) -1 



what is the mass of the cold molecular gas reservoir in these star- 
bursts? Using the IRAM Plateau de Bure Interferometer we have 
conducted a pilot study of the cold gas properties of two dusty star- 
burst galaxies in the distant cluster CI 0024+16. The initial results 
are presented here. Throughout we assume h — 0.7 in units of 
100 km s" 1 Mpc" 1 , Aa = 0.7 and f2 m = 0.3. The projected scale 
at z — 0.395 in this model is 5.3 kpc/", and the luminosity distance 
is D L ~ 2140 Mpc. 



2 OBSERVATIONS 

The two galaxies chosen for this pilot study were selected as two 
of the brightest in our Spitzer Space Telescope InfraRed Spectro- 
graph (IRS) survey of 24/im-selected members of CI 0024+16 (z = 
0.395). Full details of the mid-infrared observations can be found 
in Geach et al. (2006 & 2009). Both galaxies have S24 > 0.9 mjy, 
and implied far-infrared luminosities Lfir > lO 11 !/©; they are in 
the LIRG class. The far-infrared luminosity is estimated from the 
luminosity of the 7.7/im polycyclic aromatic hydrocarbon (PAH) 
emission line, extrapolating from the tight correlation between PAH 
strength and total infrared luminosity observed in local star form- 
ing galaxies (Smith et al. 2007; Geach et al. 2009). Both galaxies 
are in the outskirts, or infall region, of the cluster: they large clus- 
tocentric radius, 1.9-Ry and 2.8Rv, where Ry is the virial radius, 
~1.7Mpc (Treu et al. 2003). Further details of the two galaxies is 
provided in Table 1 . 

The observations were conducted using the IRAM Plateau de 
Bure Interferometer on 31 July 2008 (MIPS J002621.7) and 25 Au- 
gust 2008 (MIPS J00272 1.0) as part of program S035. The expo- 
sure times were 9.5 hr and 8.2 hr for each source respectively us- 
ing 5 antennae. The observing conditions were excellent in terms 
of atmospheric phase stability, however any anomalous and high 
phase-noise visibilities were flagged. Data were calibrated, mapped 
and analysed using the IRAM GILDAS software (Guilloteau & Lu- 
cas 2000). Phase and flux calibration was performed using 3C454.3 
and 0007+106. Secondary flux calibrators also included the sources 
MWC349, 3C454.3, 0119+115, 3C84, 2145+067. We targeted the 
CO(1-0) 1 15.27 GHz rotational transition, which at z — 0.4 is red- 
shifted into the 3 mm band, with ^ bs = 82.63 GHz. The central 
frequency of the 3 mm receiver was set to the frequency of the red- 
shifted CO(1-0) line at the systemic redshift of each galaxy de- 
rived from optical spectroscopy (Table 1). The correlator was set- 
up with 2.5 MHz spacing (2x64 channels, 320 MHz bandwidth) to 
comfortably cover the CO line should it be offset from the systemic 
redshift, and to take into account a broad velocity profile. The noise 
per 1 8 km s " 1 channel for the two observations is 1 .4 mjy beam" 1 
and 1.3 mjy beam" 1 . 



3 RESULTS & DISCUSSIONS 

3.1 CO(1-0) line luminosity and H2 gas mass 

We detect CO(1-0) emission in both galaxies with significances of 
9.3<r and 9.5<r respectively. Note that these observations represent 
the highest redshift CO detection in un-lensed LIRGs yet achieved 
(see Melchior & Combes 2008). In Figure 1 we present velocity 
integrated emission maps, and the spectra extracted from the peak 
pixel of these are shown in Figure 2. We estimate the line lumi- 
nosities by fitting the spectra with Gaussian profiles, although in 
the case of MIPS J002721.0, the broad line is better fit with a dou- 
ble Gaussian (there is a marginal improvement on the % 2 fit)- We 
discuss the physical interpretation of the line profiles in §3.2. 

We calculate L' co (in the usual units of Kkms -1 pc 2 ) from 
the integrated line emission according to e.g. Solomon & Vanden 
Bout (2005): L' co = 3.25 x W 7 Sco AvD 2 L ^ s (l + z)" 3 . The 
line luminosities are given in Table 1 with 1-cr errors estimated by 
repeatedly re-evaluating the line fits after adding noise to the spec- 
tra. The 'noise' is randomly drawn from a Gaussian distribution 
with a width equivalent to the variance in the data sampled in the 
extreme wings (|Vlsr| > 500kms _ ) of the line. We take the 
standard deviation in the resulting L' co values from 1000 of these 
'bootstrap' realisations as the la uncertainty in the measured value. 

It is assumed that the CO luminosity linearly traces the 
molecular gas mass: M(H2) = aL' co (e.g. Young & Scoville 
1991; Solomon & Vanden Bout 2005) with the conversion fac- 
tor a typically taken as either a — 4.6 Af^Kkm" 1 pc 2 )" 1 
for 'normal' spirals (the Galactic conversion factor), or a = 
0.8 Mq(K km -1 pc 2 ) -1 for ultraluminous infrared galaxies 
(ULIRGs, I/fir > W 12 Lq). The cluster starbursts are not quite 
in the ULIRG class, but are at the extreme end of what would be 
considered a 'normal' IR-bright galaxy (see Gao & Solomon 2004). 
We therefore adopt the more conservative ULIRG conversion when 
estimating the H2 gas mass (Table 1), and we discuss the implica- 
tions of this in §3.4. In this sense, the estimated H2 masses could 
be considered lower limits. 

3.2 Morphology, geometry & dynamics 

Deep, high-resolution (FWHM~0.8") optical imaging of these 
galaxies (Fig. 1, inset) allows us to constrain their orientatiorQ. 
MIPS J00262 1.7 appears to be nearly face-on to the line of sight, 
with a bright core and distinctive arms. MIPS J00272 1.0 is slightly 
more complicated; an un-sharp masked i?-band image reveals a ring 
of bright knots and a bright core or bar that appears to be slightly 

1 Unfortunately both galaxies are outside the sparse HST WFPC2 mosaic 
of this cluster (Treu et al. 2003) 
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Figure 1. Velocity integrated, clean CO(l-O) emission maps of the two galaxies in this study; each map is integrated over 210kms~ 1 . The size and shape of 
the beam is indicated by the white ellipse, and all images are orientated North top, East left. Contours are shown at integer multiples of the typical r.m.s. noise, 
~0.4mJy, in each map (negative contours down to — 2a are identified with dashed lines). The inset images are 10 X 10 arcsec cut-outs of a Subaru R-band 
image (Kodama et al. 2006) around the position of the galaxy. We present two versions of the thumbnail: the standard grey-scale luminosity (overlaid with 
CO(1-0) contours at the same levels as the main map), and an un-sharp masked version. The un-sharp masked image is made by subtracting a Gaussian kernel 
smoothed version of the R-band image (kernel FWHM corresponds to 1 kpc) from the un-processed version and reveals some of the small scale structure not 
visible in the original image. Note that MIPS J002621.7 appears to be a near face-on spiral and MIPS J002721.0 exhibits either an edge-on disc structure, or is 
made up from two components separated on scales of ~10kpc. 



extended to the north-east. It could also be orientated close to face- 
on. 

The CO(1-0) spectrum for MIPS J002721.0 exhibits a double 
peaked line profile (Fig. 2), and we show in Figure 3 that the peak of 
the CO emission traced by this velocity shear describes a locus that 
runs ~ 10 kpc across the direction of the extended optical emission. 
This could be interpreted as a rotating disc, with an angular separa- 
tion between the velocity peaks of (1.8 ± 0.8)" or (9.5 ± 4.2) kpc. 
Note that the peak of the CO emission at the [O II] redshift is cen- 
tred on the bright central core. An alternative situation is that this 
galaxy is in the process of a merger along the line-of-sight, and 
the CO(1-0) emission profile originates from emission from two 
gas-rich clumps. 

Finally, we can use the CO line widths to estimate the galax- 
ies' dynamical masses. For a rotating disc of radius R, the dynami- 
cal mass enclosed within R is Md yn sin 2 i = 4 x 10 4 i?l / p WHM 
(Neri et al. 2003). One way of estimating the likely size of the 
molecular disc is to force these galaxies to fall on the global 
Schmidt-Kennicutt relation, connecting star formation surface den- 
sity to the gas surface density: E SF oc £^ (Schmidt 1959; 
Kennicutt 1998). Assuming the gas masses and SFRs given in 
Table 1 , this implies discs of radius ~0.4 kpc and ~ 1 kpc for 
MIPS J00262 1.7 and MIPS J00272 1.0 respectively. These values 
are not unreasonable, given that observations of local LIRGs in- 
dicate large concentrations of gas in their inner few kpc (Sanders 
& Mirabel 1996). However, clearly the latter value is at odds with 



the ~ 10 kpc separation of the velocity peaks in MIPS J00272 1.0, 
and of course, these radii would further increase if we had adopted 
a higher CO-to-Fb mass conversion. Assuming i = 10° for 
MIPS J00262 1.7, we find M dy „ = (1.1 ± 0.2) x 1O 1O M . 
This ignores the systematic uncertainty on the inclination angle, 
which would clearly affect these mass estimates. The orientation 
of MIPS J00272 1.0 is similarly uncertain - it could be a face-on 
merger, or an inclined disc. Assuming it is a disc of radius 5 kpc, 
with i — 30° (appropriate for a random orientation), this yields 
M dyn — (2.9 ± 0.9) x 1O 1O M0. In a merger model this estimate 
would increase by a factor ~2 (Neri et al. 2003; Genzel et al. 2003). 



3.3 Star formation efficiency & timescales 

If we assume that the far-infrared luminosity is dominated by dusty 
star-formation, then the star formation efficiency can be represented 
by Lfir/^co- F° r these galaxies, we have Lfir/L' co = 46 ± 2 
and Lpir/Z/q = 28 ± 2. This is in the range observed in local 
(z < 0.1) spirals, which have Lfir/ 'L' co < 100 (Solomon & Van- 
den Bout 2005). In contrast, much larger values of Lfir/ 'L' co > 
100 are observed in ULIRGs and SMGs over a wide range of red- 
shift (0.1 < z < 3) implying higher SFEs, perhaps driven by more 
violent mergers and interactions (Downes, Solomon, & Radford 
1993; Tacconi et al. 2008; Coppin et al. 2008). 

Recently, Daddi et al. (2008) have presented evidence for sim- 
ilar 'low-efficiency' star formation in two BzA'-selected galaxies at 
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Figure 2. Millimeter spectra of the galaxies. The redshift of the target 
galaxy derived from their [O II]A3727 emission line is indicated by the 
vertical dashed lines, and expected to be accurate to ~30kms , The sys- 
temic CO redshifts agree with the optical redshifts (Czoske et al. 2001; 
Moran et al. 2006), although in the case of MIPS J002721.0 the broad blue- 
shifted profile is best fit by a double Gaussian with a peak separation of 
(180 ± 50)kms -1 . This could be evidence of a rotating disc, or could 
also be interpreted as a merger. Note that redshifts have been corrected for 
heliocentric velocity (see Table 1). 

z ~ 1.5. Those authors argue that while the most luminous galax- 
ies such as ULIRGs and SMGs are indeed special cases of intense, 
high-efficiency star-formation, the general case could be a less ef- 
ficient (yet still vigorous) mode of star-formation fuelled by a large 
reservoir of gas. Unfortunately a lack of CO observations of LIRGs 
at intermediate redshift makes such comparisons challenging. More 
extensive surveys in the ALMA era will no-doubt dramatically im- 
prove our understanding of this issue. 

Given our measured gas mass, we can estimate the gas deple- 
tion timescale from the galaxies' current SFRs (Table 1). Assum- 
ing the activity continues at the observed rate, and that all the gas 
is converted to stars, the lifetime of the starbursts are ~100Myr 
and ~160Myr for MIPS J002621.7 and MIPS J002721.0 respec- 
tively. Note that this is much shorter than the time it takes for a 
galaxy on a radial orbit to traverse the distance from the outskirts 
of the cluster to the core, which is of order 2-5 Gyr for CI 0024+16 
(Treu et al. 2003). The implication of these short time-scales is that 
self-exhaustion of the gas-reservoirs dominates these galaxies' star 
formation histories. Moreover it is only the environment within a 
few 100 kpc of the cluster core that can ram-pressure strip the disc 
gas (and thus truncate star formation). The major episode of star- 
formation for these very active galaxies will be long-since over by 
the time the galaxies reach the core. Note that this may not be true 
of the 'general' star forming population, where lower SFRs mean 
that gas is not exhausted in such a dramatic way, and allows star- 
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Figure 3. Un-sharp masked 10" X 10" R-band images of MIPS J00272 1 .0. 
The main morphological details are a ring of knots and a central bright 
core which is extended to the north-east. The four images are overlaid with 
the CO(1-0) emission slices taken at 70kms~ 1 intervals. The change in 
peak position represents the velocity shear across the galaxy. The contours 
are taken from the cleaned map, and are at integer multiples of the noise, 
~0.65 mJy beam - . The separation of the velocity peaks correspond to 
(9.5 ±4.2) kpc. 

formation in less active systems to endure further into the cluster 
environment (e.g. Moran et al. 2007). 

3.4 Ordinary, extraordinary or somewhere in-between? 

We have treated these galaxies as ULIRGs in terms of H2 mass 
estimation, but how would our conclusions change if we applied 
the Galactic conversion of a — 4.6? The obvious implication is 
that our derived gas masses would increase by a factor 6x, and 
so the galaxies' descendants could have stellar masses of order 
~10 11 Mq. If the majority of this additional mass is built in the 
bulge region, then this further supports a model where spiral galax- 
ies in distant clusters can evolve into some of the most massive 
SO galaxies in the cores of local clusters. The longer gas-depletion 
timescales (~500Myr and ~900 Myr) are still shorter than the time 
it takes the galaxies to reach the virial radius in CI 0024+16; our 
conclusion that the majority of the bulge growth occurs in the out- 
skirts of the cluster still holds. 

These galaxies are clearly different in nature to the Milky Way 
- they are forming stars at ~15x the Galactic rate for example. 
Although these galaxies are not as extreme as ULIRGs, it is likely 
that their mode of star formation is more similar to ULIRGs than 
the 'quiescent' mode of star formation that occurs in the discs of 
local spirals. It could be driven by a circumnuclear mode where 
cold gas is being funnelled into the central regions for example. 
Taking this into account, we assume that the ULIRG conversion 
factor is most appropriate here. 

3.5 The cluster environment and the starbursts' fate 

This pilot study is not extensive enough to make strong conclusions 
about the potential role of the cluster environment on the gas prop- 
erties of these galaxies; a more complete survey is needed. However 
we can make some comments about the possible mechanisms that 
might be operating, the likely fate of these galaxies, and the impli- 
cations this has for models of galaxy evolution in rich clusters. 

At their large clustocentric radii, these galaxies are not af- 
fected by ram-pressure stripping (although weak pressure of the 
ICM on the discs of the galaxies could promote star formation in 
the disc; see Bekki & Couch 2003 and the discussion in Geach et al. 
2009). Loosely bound gas in the halos of the galaxies could be grad- 
ually removed (Moran et al. 2007), but the additional mass such a 
reservoir could supply in the short-term is probably insignificant 
compared to the large masses already in the discs and cores. So, 
provided the molecular reservoirs are securely bound to the galaxy, 
the cold gas mass we measure now is likely to contribute to their 
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final stellar masses. Their current stellar masses (estimated from 
their if -band luminosities, see Geach et al. 2009) are likely to be 
in the range M* ~ 0.5-1.5 x 10 10 M©, implying that the cur- 
rent episode of star formation could double their stellar mass (if the 
bursts are half-way through their lifetime). Assuming most of the 
current star formation occurs in the central regions, then relative 
fading of the disc component compared to the maturing bulge will 
contribute to the morphological transformation necessary for them 
to evolve into a massive SO (Kodama & Smail 2001). 

These galaxies appear to be undergoing significant evolution 
in the outskirts of the cluster, with their activity potentially driven 
by mergers and interactions. So how important is the central clus- 
ter environment in transforming spirals to SOs, and maintaining 
the evolution of the morphology-density relationship? Our obser- 
vations suggest that the cluster environment is not required to ter- 
minate star formation in the progenitors of SOs; the key environ- 
mental effect might be more subtle. One can speculate on the likely 
processes that the galaxies will experience as they are virialised 
(long after they have used up their gas). First, the lack of accre- 
tion of additional gas from a cooling halo will prevent additional 
star formation. It is unlikely that the galaxies will experience fur- 
ther mergers, and so they cannot acquire cold gas that way; all that 
remains is morphological transformation. If the galaxies are on ra- 
dial orbits, then high speed passes in the cluster potential and close 
encounters could promote tidal disruptions or even removal of any 
remaining spiral disc (Moore, Lake & Katz 1998; Gnedin 2003), 
completing the transformation to SO. 



4 SUMMARY 

We have presented new IRAM CO detections of two LIRGs in 
the outskirts of the rich cluster CI 0024+ 16 at z = 0.395. From 
the CO(1-0) line luminosities, we measure H2 gas masses of 
(5.4 ± 0.5) x 1O 9 M and (9.1 ± 0.9) x 10 9 M Q in the two 
galaxies. This assumes a conversion between CO luminosity and 
H2 mass that is applied in ULIRGs; if we adopt the Galactic con- 
version factor, these estimates would increase by a factor ~6 x . The 
infrared-derived SFRs are ~60 M Q yr _1 , and so they will exhaust 
their reservoirs in as little as ~150Myr. 

Our observations hint that dusty starbursts in the outskirts of 
CI 0024+ 16 can assemble M* ~ 10 10 M Q of stellar mass long be- 
fore they reach the cluster core, implying that the cluster environ- 
ment is unlikely to have a significant influence on the burst. How- 
ever, it is interesting to note that one of the galaxies shows evidence 
of dynamic disturbance, potentially linked to a merger. If the galax- 
ies are destined to evolve into SOs, then the main transformation re- 
quired after the gas reservoirs have been exhausted is fading of the 
disc relative to the bulge. Potentially this could be accompanied by 
some morphological disturbance when the galaxies reach the high 
density core. 

On a final note, we highlight the fact that if the critical epoch 
of star formation occurs ~l-2Gyr before the galaxies reach the 
virial radius, and these starbursts are feeding the local SO popu- 
lation, then there should be a sizable population of post-starburst 
(k+a) galaxies at ~i? v ir in rich clusters at z ~ 0.2-0.3 (Couch 
& Sharpies 1987). Detection of such an 'intermediate' population 
would be a key piece of evidence in the model where starbursts 
triggered early during cluster infall at high-redshift are rapidly ex- 
hausting their gas reservoirs and sinking to the bottom of the cluster 
potential well where they are destined to complete their evolution 
into local massive SOs. 
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